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ABSTRACT 


Heavy quark decays provide a very advantageous investigation to test the Standard 
Model (SM). Recently, promising experiments with b quark, as well as the analysis of the 
huge data sets produced at the B factories, have led to an increasing study and sensitive 
measurements of relative b quark decays. In this thesis, I calculate various observables in 
the semi-leptonic decay process A& —y A c ru T both in the SM and in the presence of New 
Physics (NP) operators with different Lorentz structures. The results are relevant for the 
coming measurement of this semi-leptonic decay at LHC b experiment in CERN, and also 
provide theoretical predictions to refine the physics beyond the SM. 
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CHAPTER 1 


INTRODUCTION 


For the past several decades, the Standard Model (SM) has been the most successful 
theory concerning the fundamental particles and most of their interactions, namely the 
electromagnetic, weak and strong forces. It has not only successfully explained almost all 
the elementary particle experiment results so far, but precisely predicted a very wide variety 
of phenomena, leading us to a better understanding of the fundamental structure of matter. 
Specially, the discovery of the Higgs boson Bca. which is a scalar particle, makes the SM a 
remarkably successful description of the subatomic world. 

However, there are some things the SM still cannot explain, i.e., the mass of neutrino, 
the dark matter and the dark energy, and even the most familiar force in our everyday life- 
gravity. Therefore, finding physics beyond the SM becomes a major part of particle physics. 
In this explorative process of New Physics (NP), both the third generation charged leptons 
and the third generation quarks play important roles since they are comparatively heavier 
and also relatively more sensitive to NP. In addition, the constraints on NP involving the 
third generation leptons (r and v T ) and quarks (b and t) are somewhat weaker, leading to 
possible larger NP effects. 

Heavy quark decays provide a very advantageous investigation to test the SM. Re¬ 
cently, the BaBar Collaboration has reported their measurements of the ratio of the branch¬ 
ing fractions of B —y D^t~u t to B —y D^£~D(, [3j Hj: 

s 0.440 ± 0.058 ± 0.042, 

R(D*) = D J n ^4 = 0.332 ± 0.024 ± 0.018, (1.1) 

V ’ B(B -> D*+£~v t ) ’ v ; 

where l = e, fi. However, the SM predictions for R(D) and R(D*) are j3j Ej 6] 

R(D) = 0.297 ±0.017, 

R(D*) = 0.252 ±0.003, (1.2) 


which deviate from the BaBar measurements by 2cr and 2.7cr, respectively. (The BaBar Col¬ 
laboration itself reported a 3.4a deviation from SM when the two measurements of Eq. (1.1) 
are taken together.) These non-universality deviations could be providing a hint of NP 
[3 El El ITCTj . Another possible test of such a non-universality can be in the semi-leptonic 
A b —y A c tv t decay, which has not been measured experimentally though it might be mea¬ 
sured at LHC b experiment in CERN soon. In both B meson and A h baryon decays, the 
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underlying quark level transition b —>• ct~D t can be probed, as both B meson and A& baryon 
contain a b quark which will decay here. 

In this thesis, I calculate various observables in the semi-leptonic decay process A & —> 
A c tu t both in the SM and in the presence of NP operators with different Lorentz structures 
by using constraints on the NP couplings obtained by using Eq. (1.1). Since the calculations 
involve the structures of both B meson and A& baryon, the Quantum Chromodynamics 
(QCD) for the strong interactions between quarks and gluons (specially, the form factors), 
will be briefly introduced as well as the SM and the weak interactions. 


1.1 Standard Model 


The Standard Model of particle physics, formulated in the 1970s, is a theory of 
fundamental particles and their interactions. It is based on the quantum theory of fields and 
provides the most accurate description of nature at the subatomic level so far. According to 
this model, all matter is built from a small number of fundamental spin-| particles, called 
fermions', six quarks and six leptons , which follow the Fermi-Dirac statistics; while the 
carriers of the interactions are characterized as bosons , which possess integer spin (either 0 
or 1 ) and follow the Bose-Einstein statistics. There are seventeen named particles in the 
SM, which are organized in Fig. 1.1 The Higgs boson, as the last particle in the SM, was 
discovered in 2012 B 0 . 

There are four known fundamental interactions in the universe: the gravitational, 
the electromagnetic, the weak and the strong interactions. They work over different ranges 
and have different strengths. Gravity, acting between all types of particle, is the weakest 
but it has an infinite range. It is supposedly mediated by exchange of a spin-2 boson, the 
graviton, which has not been observed. Even though it is universal and is dominant on 
the scale of the universe, gravity is not included in the SM because it is much weaker than 
the other forces and can be neglected at the level of individual subatomic particles. The 
electromagnetic interaction acts between all charged particles and is mediated by photon ( 7 ) 
exchange. It also has infinite range but it is many times stronger than gravity. The weak and 
strong interactions are effective only over a very short range and dominate only at the level 
of subatomic particles. The weak interaction is associated with the exchange of elementary 
spin-1 bosons between quarks and/or leptons. These mediators are W ± and Z° bosons, with 
masses of order 100 times the proton mass. The strong interaction, as its name suggests, is 
the strongest of all four fundamental interactions. It is responsible for binding the quarks 
in the neutron and proton, and the neutrons and protons within nuclei. The strong force is 
mediated by spin- 1 , massless particles known as gluons, which couple to color charge, rather 
like the photons couple to electromagnetic charge. 

Fermions are fundamental matter particles in the SM. These twelve particles (six 
leptons and six quarks ) can be grouped into three generations. The lightest and most stable 
particles make up the first generation, whereas the heavier and less stable particles belong 
to the second and third generations. The leptons carry integral electric charge. The charged 
leptons are the electron, muon and tau, while the neutral leptons are the corresponding 
neutrinos. A different “flavour” of neutrino is paired with each “flavour” of charged lepton, 
as indicated by the subscript, i.e., (e, u e ), (//, z/ M ) and (r, u T ). The charged muon and tau 
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Figure 1.1: The SM of Elementary Particles (matter fermions in the first three generations, 
gauge bosons in the fourth column, and the Higgs boson in the fifth) 


are both unstable and decay spontaneously to electrons, neutrinos and other particles. The 
mean lifetime of the muon is 2.2 x 10 -6 s, that of the tau only 2.9 x 10 -13 s. Neutrinos were 
postulated by Pauli in 1930 in order to account for the energy and momentum missing in the 
process of nuclear /3-decay. They experience the weak interactions only. The quarks carry 
fractional electric charges, of +|e or — |e. The quark “flavour” is denoted by a symbol: u 
for ‘up’, d for ‘down’, s for ‘strange’, c for ‘charmed’, b for ‘bottom’ and t for ‘top’. While 
leptons exist as free particles, quarks are not found to do so. The bound states of quarks are 
called hadrons, which can be categorized into two families: baryons (made of three quarks) 
and mesons (made of one quark and one anti-quark). Each quark carries one of the three 
colors(or color charges): r, g and b. Quarks are bound together by gluons, which are also 
colored. Fig jl.l| shows that the three lepton pairs are exactly matched by the three quark 
pairs. 


1.2 Weak Interactions 

The weak interaction is mediated by three massive bosons, the charged W ± and 
the neutral Z°. The W + and W~ are anti-particles of each other, while the Z°, like the 
photon, is its own anti-particle. Depending on whether leptons and/or hadrons are involved, 
the weak interaction can be conventionally divided into three categories: (i) purely leptonic 
processes, e.g., fi~ —> e~ + 9 e + u^, (ii) semi-leptonic processes involving both hadrons and 
leptons, e.g., neutron /3-decay n —> p + e~ + P e , and (iii) purely hadronic processes, e.g., 
A —y p + 7T“. Perturbation theory is valid for weak and electromagnetic interactions. In 
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the 1960s, a theory of electroweak interactions was developed by Sheldon Glashow, Abdus 
Salam and Steven Weinberg that can unify the electromagnetic and weak interactions. 

So far, the experimental data on a wide range of leptonic and semi-leptonic processes 
are consistent with the assumption that the lepton fields enter the interaction only in the 
combinations 


Ja( X ) = 

l 

4M = 


(1.3) 


where J a (x ) and (x) are called leptonic currents, l = e,/i, r, i^i and i\) Vl are the corre¬ 


sponding quantized fields in Eq. (1.3). 
the transmission of quanta, i.e., 


We can describe the weak interaction as due to 


For example, the interaction Hamiltonian density of 
quantum electrodynamics (QED), according to the intermediate vector boson (IVB) theory 
can be given by 

Ui(x) = g w J a \x)W a (x) + g w J a (x)Wl(x), (1.4) 

where gw is a dimensionless coupling constant and the held W a (x) describes the W bosons 


in Eq. (1.4). This interaction is known as a “V-A” interaction, since the current ,J a (x) can 
be written as the difference of a vector part ( 7 ^) and an axial vector part (y^y 5 ). 


1.3 QCD 

Quantum chromodynamics (QCD) is the standard theory to describe the strong in¬ 
teractions, in which the color quantum number has been introduced as an extra degree of 
freedom. The color charge of a quark has three possible values, r, g and b, while anti-quarks 
carry anti-colors, f, g and b. The mediating bosons of the quark-quark interactions are called 
gluons, each carrying a color and an anti-color and postulated to belong to an octet of states. 

Quarks and gluons are observed indirectly, which means that the evidence of their 
existence inside hadrons exists but these particles have not been observed singly. Experi¬ 
ments to study the strong interactions are performed with hadrons, not with the quarks and 
gluons that are described by quantum held theory (QFT). To explore or determine the quark 
and gluon structure of hadrons, structure functions are introduced to give the properties of a 
certain particle interaction without including all of the underlying physics. The experimental 
technique is to measure the angular distribution of some processes and compare it to that 
from a point particle, then the structure of the hadron can be deduced from some form fac¬ 
tors (functions of the transferred momentum square). As an example, a charge distribution 
with electrons can be probed by measuring the cross section for scattering electrons: 

% = (1.5) 

where q is the transferred momentum and F(q) is the corresponding form factor. 
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CHAPTER 2 


FORMALISM 


The physics of the decay process A b —y can be described by an effective Hamilto¬ 

nian. In the presence of NP, the effective Hamiltonian for the qnark-level transition b —y cl~Vi 
can be written in the form mm 


n 


eff 


G F V cb 

V2 


C 7 M (1 - 75 )b + ^lC 7 m (1 - 75 )b + gRcry^l + 75 )b 1^(1 - 75 )v t 


g s cb + g P cry 5 b 1(1 - 75 )vi + h.c\, 


( 2 . 1 ) 


where G F = 1.1663787 x 10 _5 GeR~ 2 is the Fermi coupling constant, V cb is the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix element, gL,R,s,p are NP couplings and I use = 
i [ 7 ^, 7 „]/ 2 . In this thesis, I have assumed the neutrinos to be always left chiral and the 
NP effect is mainly for the r lepton. Here, I do not consider tensor operators in my work. 
Moreover, I do not assume any relation between b —> ul~ui and b —» cl~ui transitions and 
hence the analysis does not include constraints from B —y tu t . As it is expected, the SM 
effective Hamiltonian corresponds to 57 = g F = gs — 9p — 0. 

In Refs. |9], 10J, the authors had parametrized the NP in terms of the couplings gs, 
gp , gv = Qr + 9l and gA = 9r ~ 9l while in this thesis I have used g F and g F instead of g v 
and gA to align the analysis closer to realistic models [13]. The couplings gL,R,p contribute 
to R(D*) while gL,R,s contribute to R(D ). The NP couplings are considered one at a time 
and the constraints on these couplings are obtained from R(D^). 


2.1 Decay Process 


A 5 is a baryon of three quarks: u, d and b, while A c is a baryon with u, d and c quarks. 
The semi-leptonic decay process under consideration is 


Technically, the 
I 11 


A b(p) -»• T (pi) + Vt(P2) + A c (p 3 ), 

where p, pi, P 2 and ps are four energy-momentum vectors respectively, 
decay process Aft —* A c rP r transits the 5-quark to the c-quark, as is shown in Fig. |2.1 
the SM, the mediating boson is W~, which will subsequently decay into a r lepton and r 
anti-neutrino. With NP, instead of W~, the mediator can be a new particle, i.e., another 
new vector particle W'~ or a scalar (or Higgs) particle H~. I will consider these NP effects 
and compare the results with those from the SM. 
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Figure 2.1: A b Decay Process in the SM and with NP(some new mediating particles) 


2.2 Partial Decay Rate 

The partial decay rate of a particle of mass m into n bodies in its rest frame is given 
by 


dT = l^L|jvf| 2 ^ n (p;p ...,p n ), ( 2 . 2 ) 

2 m 

where M is the Feynman amplitude and is an element of n-body phase space given by 

n n ,3 

d$n(P]Pl,-,Pn) = $\P- ^Pi)IJ (2ir)32E- ' ( 2 ' 3 ) 

i=l i=l ' ' * 

This phase space element can be generated recursively 

d$n(p',Pi,-,Pn) = x d$ n _ j+1 (p; q,p j+1 , ...,p n )(2TT) 3 dq 2 , (2.4) 

where q 2 = (X^i=i -^*) 2 — ISi=iP*| J - This form is particularly useful in the case where a 
particle decays into another particle that subsequently decays, e.g., A b —> A C W~ —> A c t~D t . 
A useful method to achieve the integration of dT is given in Appendix [B} 


2.3 Feynman Amplitude 

The Feynman amplitude M includes all the physical processes. In fact, the key point 
to calculate the decay rate is to evaluate the \M\ 2 appearing in Eq. (2.2). 

The total Feynman amplitude here is 


Altotal — MsM + M„ + Mq + Mqq + M t 


l 9p 1 


(2.5) 


where M, 


9l,r,s,p 


are corresponding to the NP couplings gL,R,s, 


p- 


6 
















In the SM, the Feynman amplitude for this process is given by 


M sm = (2.6) 

where the leptonic and hadronic currents are 

£" = Mt(pi)t“(1 — TsKvte), 

= (Ac|e7„(l-7 5 )|A t ). (2.7) 

The hadronic current is expressed in terms of six form factors, 


(A c |c7^|A b ) = ma c (/i7 m + 

(A c |c 7 m 7 5 6 |)A 6 = « Ac (gi7A5 + ig2(r liV q v ls + g 3 q^ 5 )uA b - (2.8) 

Here q = p — P 3 is the transferred momentum and the form factors are functions of q 2 . 

When NP operators appear, we can obtain the hadronic current by considering the 
following relations: 


g M (A c |c 7 M 6 |A fe ) = q‘ x u Xc (hl,x + if2(7nvq v + hq»)u\ h , 

(A c |c7„75&|A & ) = q il u Xc {gil^ + m(y^q l/ ^ +gzq^)u Xb . (2.9) 

Using the equations of motion, we’ll finally get (the details are shown in Appendix |A|) 


(A c |cfc|A fe ) = « Ac (/i- 


(A c |c 7 5 6 |A fc ) = UK c (-gi 


m b - m c 
hb 


+ h- 


r 


m b - m c 


>A fe , 


m b + m c 


q 2 7s > 

3 s m b + mj At ' 


( 2 . 10 ) 


where rn b and m c are the masses of the b quark and c quark. 

To obtain the corresponding unpolarized decay rate, we need to average \M \ 2 over all 
initial polarization states and sum it over all final polarization states, 

1 


\MP = ^|M| ; 

spin 


G)\V cb \ 2 




( 2 . 11 ) 


where L^ stands for the leptonic part and H^ u stands for the hadronic part ( L^ and H^ u 
are tensors in the SM and with a vector NP effect). In the following part of this section, I 
will give the details of how to obtain \M \ 2 both in the SM and with NP effects. However, 


from Eq. (2.5) we know that there should be some crossing terms between the SM and NP 


effects. We are not going to consider the crossing term between two different NP effects since 


we just consider one NP coupling at a time. To get the final result form of Eq. (2.11), we 
have to consider the kinematics of the decay process. The kinematics here is considered in 
the rest frame of A b , and details are given in Appendix [B] 
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2.3.1 


SM 

In the SM, the leptonic tensor in Eq. (2.11) is 

^Z L SM = ^[MPl)Y(l-l5)Vv T (P2)][Vv T (P2h U ( 1 -l5)u T (Pl)} 


spin 


spin 


= Tr[{jfi + mi)7 M (l - 7 5 )(]/ 2 - m^y^l - 75)] 
= 8 (■ -g^pi ■ p 2 ~ ie^ pa pi p P2a + Pi P2 + P1P2 ) • 


( 2 . 12 ) 


Here, m\ is the mass of r lepton and m 2 is the mass of tau neutrino. I already treat the 
neutrino as massless, m 2 —y 0. The hadronic tensor in Eq. (2.11) in the SM is 


where 


Y. HsMpv 

spin 


Y^. HsmpHsmI 

spin 

— B)(A* — B*) 

spin 

J2(AA* - AB* - BA* + BB*), 

spin 


A = u Xc (filp + if2cr^q u + f 3 qp,)u Xb , 

B = U\ c (Sl7/*75 + 75 + gz%lb)u\ b i 

A* = ux^fc'yv - if 2 cr u sq s + f 3 q u )u Xc , 

B * = u\ b (< 717^75 + mcrusq 5 ^ - g^i^u^. 


(2.13) 


(2.14) 


2.3.2 NP Effects 

Let’s consider the NP effects for one NP coupling at a time and set the others to zero. 
For the vector NP effects, we consider the case with only 57 present and the case with only 
gn present. For the scalar/pseudoscalar NP effects, we consider the case with only gs or gp 
present. 

For the vector NP effect with only gp present, the NP coupling will appear in the 
hadronic current part in the Feynman amplitude and the leptonic current part is the same 
as that in the SM. Therefore, the leptonic tensor L 1 ^ = L pu M has the form in Eq. 
hadronic tensor with only gp present is 

H S M+g L ^ = Y. H S M+g LLl HsM+g L l 

spin 

= YjHsMpHsMl + H S MpHg L * v + H 9 l HsmI + H gL ^H gL l) 

spin 

= + g* L + gp+ \gp\ 2 )(A - B)(A* - B*). (2.15) 

spin 


(2.12). The 


With the gp present, we will have = Lg U M for the same reason as with gp. The corre- 





sponding hadronic tensor is 


HsM+qr^ - ^2 H SM+g Rfl HsM+g R t 

spin 

= YXHs Mlt Hs M l + H SMfl H git l + H gRfi H SM t + H gR ^H gRv ) 

spin 

= Y^K A - B )( A *- B *) + 9*r( A ~ b )( A * + b *) 

spin 

+g R (A + B)(A* - B *) + \g R \ 2 (A + B)(A* + B *)]. (2.16) 


The A, B , A* and B* appearing in Eq. (2.15) and Eq. (2.16) are given in Eq. (2.14). 

For the scalar and pseudoscalar NP effects, the Feynman amplitudes are given by 

GpVcb 


Mg s — 


9 s- 


\/2 


- L S C , 


vr _ GpVcb T n 

M gp ~ 9p ^ L S D i 


(2.17) 


epto nic current, C and D are the corresponding hadronic currents. Then 


where Lg is the 
— 2 

the \M\ in Eq. (|2.11|) with NP effect with only g$ present becomes 


W SM+as = \Y,(\ M SM \ 2 + M S M M ; t + M K M' S M+\M 3s \ 2 ) 

spin 

G 2 f \V cb \ 2 ^ 

= ' J2i LL *( A - B )( A * - B *) + g*s LL *s( A - B )G* 

spin 

+g s L s L*C(A * - B*) + \g s \ 2 L s L* s CC% (2.18) 

— 12 

and the \M\ with only gp present becomes 


\M\ 2 sm +sp = + + 

spin 

G 2 f \V cb \ 2 ^ 

= 1 ' - ^[LL*(H - B){A* - B*) + g* P LL* s (A - B)D* 

spin 

+g P L s L*D(A* - B*) + \g P \ 2 L s L* s DD*}. (2.19) 


Here, L is the leptonic current in Eq. (2.7) and L* is its conjugate part. In addition, 


L s = «r(Pl)(l-75K T (P2), 


C = u Xc (h- 


D = u Xc (-g 1 


m b - m c 
hb 


+ / 3 - 


q " 


rn b - m c 


~)u\ b , 


m b + rri, 


9 2 75 N 

- . 9.3 --- )U Xb , 

m b + m c 


( 2 . 20 ) 
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and their corresponding conjugate parts are 


T* 

C* 

D* 


^ t (P 2)(1 +75)«r(Pl), 


b {f 


+ h~ 


(f 


1 - 

rn b — m c mb — m c 


-)u\c 


, _ _ , g 2 75 > 

UXb[ ' 9 l m b + m c 93 m b + mJ UXc ‘ 


( 2 . 21 ) 


The cross terms between two different NP couplings are zero since we consider one 
NP coupling at a time, as Eq. (2.5) indicates. 


2.4 Observables 


The calculation is based on integration of Eq. (2.2), which gives us the decay rate of 
the process A b —> A c t9 t directly. In this thesis, we’ll define the following observables. 


BR[A b —» A c ru T \ 
BR[A b —> A c lD(] 


( 2 . 22 ) 


Here £ represents /i or e. The branching ratio BR for a specific decay process is defined by 


BRi 



(2.23) 


where r,; is the decay rate for this process and T?: is the total decay rate. 

The differential distributions with respect to the transferred momentum square q 2 will 
be shown in the results ( dY/dq 2 ). Also, we will define the ratio of differential distributions 


B K {q 2 ) 


dr[Ab —>■ A c ru T \ 
dq 2 


/ 


dT[A b —> A c £v(\ 
dq 2 


(2.24) 


The results will show that these observables are not very sensitive to variations in the 
hadronic form factors. 
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CHAPTER 3 


NUMERICAL RESULTS 


In this section, I will present the constraints on the NP couplings, then I will discuss 
and show the form factors used in this work. Finally, I will present the result graphs of the 
observables defined in section 12.41 


3.1 NP Couplings 


The NP constraints on the NP couplings are obtained from R(R)W) in Eq. (1.1) and 


Eq. (1.2). The relative decay processes are 


B —y Dt v t , and B —y D*r u T . 

The two main reasons we use the NP constraints from the decay process B —» D^t~V t here 
are: the experimental results of R(D^) from the BaBar Collaboration deviating from those 
in SM provide a hint of NP [3- 4j, and both A^ baryon decay and B meson decay involve the 


transition b —>• cl up as shown in Fig. 2.1 and in Fig. 3.1 The formalism to constrain the 
NP couplings here is 


R 


exp 


P f 1 In m npM* sm MsmM^p , > 2 \Mnp\\ (o n 

- RSMfl + gNP— TTV Pf- + 9nP~ TT7 E - + \9NP\ 7TV wj, (.3.1J 

\Msm\ \Msm\ \Msm\ 


\Mnp\ 


where gjyp stands for gL,p,s,p- 


Ef 




Figure 3.1: B —> r u T Decay Process 
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Re(9 s ) 


Re(g P ) 


Figure 3.2: The figures show the constraints on the NP couplings taken one at a time at 
the 95% CL limit. When the couplings contribute to both R(D) and R(D*) the green areas 
indicate constraints from R(D) and the orange ones from R(D*). 


The NP coupling g$ only contributes to R(D) and gp only contributes to R(D*) while 
g^p contributes to both R(D) and R(D*). These couplings can be derived from Eq.(3.1) 
and their allowed regions are shown in Fig. |3.2[ To get them, we need the form factors of 
B —» Z)(*) decay, given in the Appendix [c| and D 


3.2 Form Factors 

One of the main inputs in our calculations are the form factors. As first principle, 
lattice calculations of the form factors are not yet available. The form factors we use here are 
from QCD sum rules, which is a well known approach to compute non-perturbative effects 
like form factors for systems with both light and heavy quarks nuns]. Another approach 
to investigate the heavy quark systems involving quark flavour and spin symmetry ,called 
Heavy Quark Effective Theory (HQET) [16]. The form factors we choose are consistent with 
HQET. 

There are different values of the relative parameters for the kinematic region ascribed 
to the continuum model. This gives us some choices of the form factors. In Ref. [15] various 
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parametrizations of the form factors are given. They are shown below. 


continuum model 

hi 

F? («)=/, 

F?(t)(GeV -') = h 

rectangular 

1 

6.66/(20.27 -t) 

—0.21/(15.15 — t) 

rectangular 

2 

8.13/(22.50 -t) 

—0.22/(13.63 — t) 

triangular 

3 

13.74/(26.68 -1) 

—0.41/(18.65 — t) 

triangular 

4 

16.17/(29.12 — t) 

—0.45/(19.04 — t) 


Table 3.1: Various Choices of Form Factors (t = q 2 ) 


The form factors in Table 3T are based on four continuum models indicated by 
k — 1,2, 3,4. However, the differences of the results from these models are very small, which 
can be shown in the result graphs. Moreover, these form factors in Table 3.1 satisfy the 
HQET relation in the mb —> oo limit. They have the following relations: 


/i — 9h — 92, fz 93 0. 


(3.2) 


3.3 Result Graphs 


I have used the following masses in my calculations. The masses of the particles are 
m = 5.6195 GeV, m T = 1.77682 GeV, m^ = 0.10565837 GeV, m 3 = 2.28646 GeV, mb = 4.66 
GeV and m c = 1.275 GeV (TTj). In the following I will present the results for observables R\ b , 
dT / dq 2 and B\ b (q 2 ). For the Erst and third observables I use different models of the form 
factors given in Table 3.1 For the differential distribution dT/dq 2 , I present the average 


results over the form factors. 


continuum model 

1 

2 

3 

4 

Average 

Ref. m 

Ref. j2U] 

RaASM) 

0.31 

0.29 

0.28 

0.27 

0.29 

0.29 

0.31 


Table 3.2: Values of R^ b in the SM 


In Table 3.2 the prediction for R\ b in the SM are given for the various choices of the 
form factors in Table 3T I also compare our results with other calculations of this quantity 
by other groups using different form factors. The average value we found for R\ b in the SM 
is R\ bt sM = 0.29. This agrees very well with values for this quantity obtained in Ref. [ T9] , 
which uses a covariant confined quark model for the form factors, and Ref. [20] which uses 
the form factor model in Ref. 123 . These results indicate that the ratio R\ b is largely free 
from form factor uncertainties making it an excellent probe to find new physics. 


Now I give the discussions of the results. From Eq. (2.5), we can make some general 


observations. We start with the case where only g Jj is present. In this case the NP Feynman 
amplitude has the same structure as the SM one and the total Feynman amplitude is just 
the SM amplitude modified by the factor (1 + g L ). Hence, if only g L is present, then 


-R/\ 


— sm 11 + 9l\ 


(3.3) 


Therefore, in this case, R\ b > R\ b Q M and we find the range of R\ b to be 0.31 ^ 0.44. The 
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Only g L Present 


Only g R Present 



Model 

Only g L Present 


' -SM 




-NP: g R =0.083-l*0.729 




' -— NP:g R =-0.1111+1*0.4961 




■ 




♦ 


■ 




♦ 




* 

- 

1 2 

3 

4 



Model 

Only g R Present 




Only g L Present 


Only g R Present 




Figure 3.3: The graphs on the left-side (right-side) show the compared results between the 
standard model and new physics with only g bj (gn) present. The top and bottom row of 
graphs depict R\ b = BR[A b — >• A c tv t ]/ BR[A b — > A c iv(\ and the ratio of differential distri¬ 
butions B\ b (q 2 ) = A b —> A c tv t ) /^ (A(, — > A c ivt ) as a function of q 2 , respectively for 

The middle graphs depict the average differential 
Some representative values of the 


the various form factors in Table 3.1 
decay rate with respect to q 2 for the process A b —> A c tv. 
couplings have been chosen. 
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Only g s Present 


Only g s Present 



■ SM 

■ NP: g s =-0.6824+l*0.7952 

■ NP: g s =- 1.03-1*0.902 


Model 

Only gs Present 


2 3 

Model 

Only g s Present 



Only g s Present 


Only g s Present 




Figure 3.4: The figures show the compared results between the standard model and 

new physics with only gs present. The top and bottom row of graphs depict R\ b = 
BR[Ab —> A c Tv r \/BR[Ab — » A c lvi\ and the ratio of differential distributions B\ b (q 2 ) = 
^2 (A b —> A c tu t )/^ (A;, —> A c tv() as a function of q 2 , respectively for the various form 
factors in Table 3.1 The middle graphs depict the average differential decay rate with re¬ 
spect to q 2 for the process A b — » A c tv t . Some representative values of the couplings have 
been chosen. 
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Only g P Present 




Only g P Present 



Only g P Present 



Model 


Only g P Present 



Only g P Present 



Figure 3.5: The figures show the compared results between the standard model and 
new physics with only gp present. The top and bottom row of graphs depict R\ b = 
BR[Ab —> A c tu t ]/BR[Ab —>■ A c Zvg\ and the ratio of differential distributions B\ b (q 2 ) = 
^ 2 (A b —> A c TV T )/j^j(Ab — > AJ'Vji) as a function of q 2 , respectively for the various form 
factors in Table 3.1 The middle graphs depict the average differential decay rate with re¬ 
spect to q 2 for the process A b — » A c tv t . Some representative values of the couplings have 
been chosen. 
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shape of the differential distribution dr/ dq 2 is the same as in the SM. In the left-side figures 
of Fig. 3.3, we show the plots for R\ b , dT/dq 2 and B Rb (q 2 ) when only gp is present. We then 


consider the case where only g R is present. If only g R is present, then from Eq. (2.13) and 


Eq. (2.16) we can deduce that no clear relation between Rp b and R A b sM can be obtained in 


this case. However, for the allowed g R couplings, we find R\ b is greater than the SM value 
and is in the range 0.30 ^ 0.51. The shape of the differential distribution dT/dq 2 is also 


the same as that of the SM. In the right-side figures of Fig. 3.3, we show the plots for Rp b , 
dT/dq 2 and B\ b (q 2 ) when only g R is present. 

We now move to the case when only gs.p are present. Using Eq. (2.5), Eq. (2.6) and 


Eq. (2.17), we can write 


Ra„ — R A bS M + \9s\ 2 As + 2Re(g s )B s , 
R^bSM + I 9p\ 2 Ap + 2 Re(gp)Bp, 


R A h = 


(3.4) 


since the physics of the the decay process only underlies in the Feyman amplitudes. The 
quantities As t p and Bs,p depend on masses and form factors and they are positive. Hence, 
for Re(gp) > 0 or Re(gs) > 0, R\ b is always greater than or equal to Ra bSM - But, for 
Re(gp) < 0 or Re(gs) < 0, R\ b can be possibly less than the SM value. However, for the 
given constraints on gs here, we can make Rp b only slightly less than the SM value while 
for gp it is always larger than the SM value. We find R\ b is in the range 0.28 ^ 0.36 when 
only gs is present and in the range 0.30 ^ 0.42 when only gp is present. In Fig. 3 A we show 
the plots for R\ b , dT/dq 2 and B\ b (q 2 ) when only gs is present. The shape of the differential 
distribution dT/dq 2 can be different from that of the SM. In Fig. 3.5 we show the plots 


for R\ b , dT/dq 2 and B^ b (q 2 ) when only gp is present. In this case also the shape of the 
differential distribution dT/dq 2 can be different from that of the SM. 


NP 


Rh. b,max 

Only g L 
Only g R 
Only g s 
Only g P 

0.31, g L = -0.065 + 0.447 i 
0.30, g R = -0.033 + 0.119 i 
0.28, g s = -1.442 

0.30, g P = 0.587 

0.44, g L = -0.144 + 0.903 i 
0.51, g R = 0.182 + 0.914 i 
0.36, g s = 0.443 

0.42, g P = -5.859 


Table 3.3: Minimum and Maximum Values for the Averaged R\ b . 


In Table 3.3, we show the minimum and maximum values for the averaged R\ b with 
the corresponding NP couplings. 


17 
















CHAPTER 4 


CONCLUSION 


In this thesis, I calculated the SM and the NP predictions for the decay A& — y A c tv t . 
Motivation to study this decay comes from the recent hints of lepton flavor non-universality 
observed by the BaBar Collaboration in R(D^) = nM+z--^ = CaO- I used a 


B(B^DM+e~9c) 

general parametrization of the NP operators and fixed the new physics couplin gs fro m the 
experimental measurements of R(D ) and R(D*). The predictions for R\ b (Eq.(2.22)), 
and B\ b (q 2 ) (Eq.(2.24)) are made by taking one of the various NP couplings at a time. We 
found the interesting results that gL,n,p couplings gave predictions larger than the SM values 
for all the three observables while the g$ couplings gave predictions which could be larger or 
smaller than the SM values. 

This thesis is related to our recent work of A& —y R c rv r decay 
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APPENDIX A: NP OPERATORS EXPRESSED IN TERMS OF FORM FACTORS 
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If we consider the hadronic current: 


(A c |c 7 M 6 |A fe ) = wa c (/i7a» + if2<*pu<f + h%)u\ b . 


(A.l) 


then (q = p — p 3 ) 


g M (A c |c 7 M 6 |A 6 ) = g M nA c (/i7M + + hq»)u Xb . 


(A. 2 ) 


The Left-Hand-Side of Eq. (A.2) is 


g M (A c |c 7 M fe|A 6 ) = (A c |cg M 7 M 6 |A b ) 

= (A c |c^ 6 |A b ) 

= (A c \c(jj> - p/ 3 )b\A b ) 

= (m b -m c )(A c \cb\A b ), 


(A.3) 


where I used the equation of motion: fib = m b b and cpj 3 = m c c, while the Right-Hand-Side 
of Eq. (1A.2I) is 


9 m Ma c (/i7m + if 2 (y^q u + hq»)u Xb = u Xc (fii + 0 + fzq 2 )u Xb . 
Thus, we can get: 


+ /s- 


r 




then 


(A c |c 6 |A b ) = u Xc (fi -A , JO 

mn b — m c m b — m c 

Now, consider 

(A c 1 07 ^ 756 ) Aft) = 'u Ac (fi'i7 M 7 5 + ig2(7^q u ib + g?,q^)u Xbl 
(A c | C7 M 7 5 6| Aft) = g^UAc (^17^75 + + g^lb)u Xb . 


The Left-Hand-Side of Eq. (A.7) is 

(A c 1 07^756] Aft) = (A c |cg /i 7 M 7 5 6|Aft) 

= (A c |c^ 7 5 6 |Aft) 

= (A c |c(^ — ^3)756 1 Aft) 

= -(m b + m c )(A c \c'y 5 b\A b ), 

Where I used the equation of motion: fib = m b b and cp/ 3 = m c c, and ^75 = 
Right-Hand-Side of Eq. (A.7) is 

^A c (5h7//y5 + mcr^q 1 '^ + g^lb)u Xb = u Xc (gip b + 0 + g 3 q 2 ^)u Xb . 
Thus, we can get: 

(A c |c 7 5 6 |Aft) = u Xc (-g 1 —^ - g 3 )u Xb . 

m b + m c m b + m c 


(A.4) 

(A.5) 

(A. 6 ) 

(A.7) 


(A. 8 ) 

The 

(A.9) 

(A. 10 ) 
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APPENDIX B: KINEMATICS 
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In the rest frame of A&, we have: 


P = (m, 0,0,0), 

Pi = (Ei,pi), 

P 2 = ( E 2 ,p 2 ), 

P3 = (E 3 ,p 3 ). (B.l) 

The transferred momentum q = pi + p 2 = p — p 3 . We have p 3 = p — Pi — p 2 - By 
considering the Lorentz invariance, we can find out the following kinematic relations: 


p 2 

= m 2 , 

2 

2 

Pi 

= m,, 

pi 

= o, 

2 

2 

P3 

= m 3 » 

p-Pl 

= mEi, 

P'P2 

= mE 2 , 

CO 

= mE 3 , 

p-q 

= mf?i + m^2 

Pi ■ q 

= ^ 2 +™?)> 

Pl ' p2 


CO 

= mEi — -m 2 

P 2 ■ q 

= ^(? 2 -mi), 

to 

CO 

= ml?2 + ^m 2 

P3 ■ q 

= ^(9 2 -«r3). 


(B-2) 


To achieve the integration of differential decay rate, let’s define p^ = pi + Pj and 
m ij = P%- Then rn 2 i 2 + m 2 23 +'m 2 13 = m 2 +m?+m 2 +m| and m 2 2 = (p-p 3 ) 2 = m 2 +m 2 3 -2mE 3 , 
where E 3 is the energy of particle 3 in the rest frame of m. 

From m 23 = (p — Pi) 2 = rn 2 + m\ — 2mE \, we have 


E\ 


m 2 + rn 2 — m 23 
2m 


(B.3) 


From m \ 3 = (p — p 2 ) 2 = m 2 + m 2 — 2 m E 2 and m 2 2 + m 23 + m \ 3 = m 2 + m\ + rn 2 + m 2 , 
we have 


E 2 


rn\ 2 + m 23 — m\ — m 2 
2 m 


(B.4) 
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From m\ 2 = (p — P 3) 2 = iri 2 + m\ — 2mE 3 , we have 


E 3 


m 2 + m 3 — m i 2 

2 m 


Using the standard form for the Dalitz plot, we can get 


dr = ~^—\M\ 2 d& n (p]pi, ...,p n ) 

2m 

= jk ^ mdEliE2 


(B.5) 


(B.6) 


Here, the Dalitz plot: for a given value of m 2 l2 , the range of m 23 is determined by its values 
when p 2 is parallel or anti-parallel to p 3 : 


(m 2 23 ) mx = (El + Eff - ( yjEf-ml - ^Ef - m|) 2 , 

04 )™.. = (K + Elf - ( s]e f - m\ + s]Ef - ml f, (B.7) 

where E 2 = {jn\ 2 ~ m \ + r?7|)/2?77i2 and E^ = ( 777 2 — rn\ 2 — m^)/2m 3 2 are the energies of 

particles 2 and 3 in the 77712 rest frame. Since m\ 2 = q 2 , the differential decay rate with 

respect to dq 2 is: 

W = (R8) 
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APPENDIX C : B ^ D* FORM FACTORS 



C.l B —» D*r Vf, Angular Distribution 


The full B —>■ D*t v tJ angular distribution is given by [9] 
r/T D * r 

= N\p D *\ 2\'H 0 \ 2 sin 2 Q x + (\H\\ | 2 + |'Hj_| 2 )(l + cos 6 Q 2 — 4i?e[*4.||'H^] cos^ 


dq 2 d cos 9[ 


2 

+-wy ^ 2 |"Ho cos 6 ^ — TLtp | 2 + (|'H ||| 2 + \R±\ 2 ) sin 2 O^J , 

where 9i is the angle between the D* meson and the r lepton three-moment a in 
frame, N = ^1 — -jr'j and the amplitude 'Ho,\\,x,t,p are given in Sec. 

definition of T-i t p is 


C.4 


U t P = H 




n f 


m-r 


(C.l) 


he q 2 rest 
Also, the 


(C.2) 


C.2 B —>• D*t iS/, Form Factors 


The relevant form factors for the B —>• D* matrix elements of the vector = c'y^b 
and axial-vector A M = 07^756 currents are defined as 


{d’K\b) = vrfw'M. 

% + Trijj* 

(D^A^B) = 2in D *A 0 (q 2 ) ( —^-q fl + {m B + m D *)Ai(q 2 ) 


e ■ q 
q z 


-A 2 {q 2 ) 


e • q 


{m B + m D * 


{jpB + PD*)p ~ 


m B — m 


D* 


% 


(C.3) 


In the Heavy Quark Effective Theory (HQET), the form factors in Eq. (C.3) are given 


by 


A 0 (q 2 ) 

^oH, , x 

= rj ^(™), 

Kd* 


Ai(q 2 ) 

- Rd* 2 h Al (w), 


A 2 (q 2 ) 

r 2 (w) 

= JD ^(™)> 

Kp* 


V{q 2 ) 

R iM, f x 

= p 

Kp* 

(C.4) 


where Rp>* = 2 ^JnnBm * D /(rag + m * D ). The summary results of w dependence of the form 
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factors mm are 


h Al (w) 

Ri(w) 

R 2 (w) 

R 0 {w) 


h Al ( 1) 1 - 8 p 2 z + (53 p 2 - 15 )z 2 - (231 p 2 


i?i(l) - 0.12(iy - 1) + 0.05(w - l) 2 , 
R 2 ( 1) + 0 . 11 (w - 1 ) - 0.06(w - l) 2 , 
i? 0 (l) - 0 . 11 (w - 1 ) + 0 . 01 (w - l) 2 , 



(C.5) 


where z = (y/w + 1 — y/2)/{y/w ± 1 + y/2). The numerical values of the free parameters p 2 , 
h Al ( 1 ), i?i(l) and i? 2 (l) are (25| 


h Al (l)\V cb \ 

P 2 

fli(l) 

H 2 (l) 


(34.6 ± 0.2 ± 1.0) x 10“ 3 , 
1.214 ±0.034 ±0.009, 
1.401 ±0.034 ±0.018, 
0.864 ±0.024 ±0.008, 


(C. 6 ) 


and i?o(l) = 1-14 is taken from Ref. |5jj. In the numerical analysis, we may allow 10% 
uncertainties in the i?o(l) value to account higher order corrections. 


Therefore, in the HQET the amplitudes in Eq. (C.l) become 


, . , rjw 2 — 1) r . m 2 R (l ± r 2 — 2r*w) 

Rtp = ~*m B (l + , 2 ± . J i Al (w)R. 0 (w)[(\ -g A )~— ^ ^ V ], 


H 0 = 


1 ± r 2 — 2 r*w 
m B ( 1 - r*)(iy ± 1 ) A /fy 


\J (1 ± r 2 — 2 r*w) 


h Al (w) 


1 ± 


(w - 1)(1 - H 2 (w)) 


(1 


mi(m b ± m c ) 
(1 - Pa), 


H|| = m B \/2n(w+ l)h Al (w)(l -g A ), 


Ui_ = m B y/2r*(w 2 - l)h Al (w)R 1 (w)(l ± g v ), 
where r* = m B */m B . 


(C.7) 


C.3 B —>■ D* Kinematics 

Here, for a = (ao, a±, a 2 , 03 ) and b = (bo, 61 , b 2 , 63 ), we have a-b — aob 0 — (ci\bi + a 2 b 2 ± 

0363). 

In the rest frame of B: 


Pb = (m B , 0,0,0), 

Pd * = (-Ed*, 0, 0, \p B * |), 

q = (g 0 , 0,0, —|pd* |), (C. 8 ) 
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where q = Pb = Pd* and 


E d » = {m 2 B + m 2 D * — q 2 )/(2m B ), 
l„- I - . l, m B + m D- .,2 


m D*i 


% = ( m% - m 2 D , + q 2 )/(2m B ). 


(C.9) 


The polarization vectors of D* are given by 


eo — - (\p D *\,0,0,E D *), 

171 £>* 

e± = T^(0,l, ± b°). 

The polarization vector of virtual gauge boson are given by 


(C.10) 


= -y=(|PD*|,0,0, -g 0 ), 

e± = ^(0,±l,-i,0), 

£t = -^=(q o ,0,0,-\p D *\). 

To get the amplitudes, we have the following useful relations: 
e* 0 -q = 0 , 

1 1 

e* 0 -p D * = ~^(\pd*\E d * + \p D *\qo) = -J=rn B \p D *|* 


(C.ll) 


(C.12) 


* 1 n IT? I I | \ m B\PD*\ 

e 0 -q = -( \pd*\E d * + \p D * g 0 ) = -, 

m D * m D * 

1 „ „ „ 


f* . <r* — 

tn tn — 


m D * \f([ 1 


(\pd*\ + E D *q 0 ), 


(C.13) 


h • ( i = 

e* ■ £* = - (\pd*\E d * + \p D *\qo) = - - m B \p D *\, 


(C.14) 
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el 



e_ 


• e_ 


^( 0 , 1 , 10 ) • 
^(0,-1, *,0) 


~(0,l,i,0) = 


= -i(0-l-l-0) 

^(0 + 1 + 1 + 0 ) = 1 , 


1, 


e* ± -q = 0. 


(C.15) 

(C.16) 


C.4 B —>■ D* Amplitudes 


In Ref. [ID], Eq. (A.6) (the V-part) used e^ vpa . I used e 0 i 23 = 1, which I think should 
agree with in Ref. [23] where the authors used e 0123 = — 1. 


1. ftn 


Vo = 0, 

A 0 = £* 0 (D*\A p \B) 

= (m B + m D *)A 1 (q 2 )e* 0 ■ £* 0 - 


A '<? 


ms + rn D * 


A 2 (q 2 )}(2p D * -£* 0 ) 


4m 2 B \p L 


~~ -t=[(™| - m D* ~ Q 2 )i.m B + m D *)Ai(q 2 ) - AAAAAA A 2 (q\$ r 17) 

2 m D * yjq 2 m B + m D * 


B 0 = (V 0 -A 0 )(l-g A ) 

:[{m% - m 2 D * - q 2 )(m B + m D *)Ai(q 2 ) 


2 m D . a Jq 2 
Am 2 R \p 12 


'D* 


m B + m D * 


A 2 (q 2 )\(l — g A ). 


(CAS) 


Here, \p D *\ 2 + qoE D * = 

2 . Up 


9 9 9 

_ m z B —m z D * —q z 


V t = 0, 

At = £* t (D*\A,\B) 

2 m D *A 0 (q 2 )e* t ■ q 


q 

Jl \ * 


m’Kfi'fi - e> q ' 


v 7 ?' 


A 2 (q 2 )e* ■ q m 2 B — m 2 D * mr B — m 2 D « 


m B + m D * 


Vq 2 


\fA 


m B \p D *\A 0 (q 2 ), 


(C.19) 
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2 


(C.20) 


/ = (V t -A t )(l-g A ) = 




m B \p D *\A 0 (q 2 )(l — g A ). 


vv = 


? + (D*\V,\B) 
2 iV(q 2 ) 
m B + m fr 
2 iV(q 2 ) 


—*u p 8 

£^p&£+e Pd*Pb 


m B + rrirr 
2 iV(q 2 ) 


(^ 1230 ^+ + ^ 2130 ^+ e * )Pd*Pb 

[—T^(-T^) + -4(“i)]|Pr»*|mfl 


*4 L = 


m B + m D * 1 ^2 y/2 V2 

2\/(g 2 )m B |p g * 
m B + m D * 
t + (D*\A,\B) 


= (m B + m jD *)44i(g 2 )e+ • £+ 

= (m B + mD*)44i(g 2 ), 


4s(g )e* • g 

m B + mg?* 


* 0 


— V + (l + g v ) — A + (l — g A ) 
2V(q 2 )m B \p D 


rn B + m D * 


(1 + g v ) - (m B + rriD^A^q )(1 - g A ). 


V _ = t_{D*\V^\B) 
2iV (g 2 ) 
mg + m D . 

2iV(q 2 ) 
mg + m D . 


— LL *V p 8 

£,wps£- e p p D *p B 


{£i23o£*-£*~ + £2i3o£t 2 £* i )p B *P l B 


—*2 *1\ 3 „0 


2 iV(q 2 ) , 1 i i 1 , 




m B + m D A y/2 ^2 V2 

2\/(g 2 )mB|p D 

m B + m d* 

= e*_(D*\A,\B) 


= (m B + mD*)44i(g 2 )el • el 
= (m B + mD*)44i(g 2 ), 


2 \ * ^ A 2 (q 2 )e*-q 


m B + mg)* 


* 0 


y_(l + 5»y) - A-{ 1 -5-a) 

2\/(g 2 )mgj|pg) 


(1 + 9V) - (m B + m D *)Ai(g L )(l - 5 -a)- 


(C.21) 

(C.22) 


(C.23) 

(C.24) 


m B + m D * 



Here, £1230 — — 1 , £2130 — 1 - 


4. n 


(ll.-L) : 


n± = 




n. 


m B + mp>* 


(C.25) 


5. H P : 


11 v 7 ^ 

= + m D *)Ai(q 2 )(l - g A ). 


q»(D*\c^ l5 b\B) = (D*\c^ 5 b\B) 

= {D*\c(p/b ~pfi*)i 5 b\B) 

= ~(m B + m D *)(D*\c'i 5 b\B). 

Here, I already used the equation of motion. Also, 

q tl {D*\c 7 ll ^ b b\B) = 2m D *\p D *\A 0 (q 2 )e* ■ q 
= 2m D *\p D *\A 0 (q 2 ) 


2 \ m B \pd* 


(C.26) 


(C.27) 


m D * 


(C.28) 


H P = {.D*\c lb b\B)g P 
-2m B \p D *\A 0 (q 


m B + rn D * 


~9p- 


(C.29) 
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APPENDIX D: B ->• DtP t FORM FACTORS 
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D.l B —> Dtv t Angular Distribution 


The B —>■ Dtu t angular distribution (the differential decay rate) for the lepton hclicity 
A T = ±4 are 


dr D [A T = -1/2] 

dq 2 d cos 61 
dr D [A r = 1/2] 
dq 2 d cos 9 1 


2 N\p D \\H 0 \ 2 sin 2 9 h 

2 

TTi 

2N\p D \^\H 0 cos9 l -H tS 

q 2 


(D.l) 


The differential decay rate corresponding to the hclicity A T = 1/2 vanishes for the 
light leptons (e,/i) as their mass is much smaller. 


D.2 B —> Dtu t Amplitudes 

The methods to get the amplitudes from the B —>• D matrix elements is the same as 
for B —>■ D*. Here, I didn’t include the details. 

The amplitudes are 


H 0 = 
H t = 
H s = 


2 m B | p 


'D 


V7 


F +(« 2 )(l + 9v). 


m r — m 


D 


v7- 


m r — m 


D 


m h - m c 


Fo(q 2 )(l + gv)-, 


F 0 (q 2 )gs- 


(D.2) 


Also, the dehnition of H t s is 


H tS = tf t + — 


m T 


(D.3) 


D.3 B —> Dtv t Form Factors 


The form factors F + (q 2 ) and F 0 (q 2 ) of the B —> D matrix elements are defined as 


(D\cYb\B) 

(D\cb\B) 


= F + (q 2 


Pb+Pd 


m n — m 


D q » 


+ F 0 (q 


m | - m- D ^ 


rnh — m 


D 


m b (p) - m c {p) 


Fo(q 2 


(D.4) 
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In the heavy quark effective theory, the form factors in Eq.(D.4) are given by 


F + {q 2 ) 

Fo(q 2 ) 


Vi(w) 

Rd ’ 

(1 + w)R D 

-O- S l( w ) 


(D.5) 


where Rd = ‘l^/ msvnD f {jns ± wd) and r = mo /tub- The parametrization of the form 
factor Vi (w) is given by 123 

Vi(w) = Vi(l)[l — 8 pfz ± (51Pi — 10)z 2 — (252p 2 — 84) z 3 ], (D.6) 

where z = (y/w ± 1 — a/2) / (\Au ± 1 + y/ 2 ). The numerical values of the free parameters are 

m 


yi(l)|Vrf,| = (43.0 ± 1.9 ± 1.4) x 10“ 3 , 
p\ = 1.20 ±0.09 ±0.04. 


(D.7) 


The parametrization of form factor S\(w) is given by [6j 

S^w) = 1.0036[1 — 0.0068(u> — 1) ± 0.0017(tc — l) 2 — 0.0013(tc — l) 3 ]Vi(tc). (D.8) 


In the HQET, the amplitudes in Eq. (D.2) becomes 


Ho 

H ts 


r(w 2 — 1) 


m ° ( - 1 + r ^(l+r*- 2 rw) 
mg( 1 — r)y/r{w ± 1) 


Vi[w](l +g v ), 


\J (1 ± r 2 — 2 rw) 


S'iH (1 + gv) ± 


m 2 B { 1 ± r 2 — 2 rw) 
mi(m b (p) - m c (/i)) 


9s 


(D.9) 
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